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SUMMARY 



Two topics of interest. to persons attempting to apply 
the heat method of preventing ice formation on aircraft 
are considered. Surfaces moving through air at high speed 
are shown, both theoretically and exper imen t ally , to he 
subject to important aerodynamic heating effects that will 
materially reduce the heat required to prevent ice. 

Numerical calculations of the paths of water drops 
in an air stream around a circular cylinder are given. 
From these calculations, information is obtained on the 
percentage of the swept area cleared of drops. 

INTRODUCTION 



The NACA ha.s been conducting a program to develop a 
method of using engine-exhaust heat to prevent or to re- 
move ice formations on aircraft. The chief object of this 
paper is to provide information that will be useful to 
persons attempting to apply this method. 

Two subjects are included in the present investiga- 
tion: 

1. The temperature rise of surfaces exposed to 

boundary-layer friction has been studied. 
This natural temperature rise is expected 
materially to reduce the amount of heat re- 
quired to keep a surface in high-speed flight 
above the freezing temperature. 

2. The paths of water drops in an air stream dis- 

turbed by the presence of a circular cylinder 
have been numerically calculated. The calcu- 
lation of these paths leads to information 
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on the percentage of the' swept area cleared 
of drops "by the. circular cylinder. This in- 
formation should "be of assistance in estimat- 
ing, the- heat required to overcome the super- 
cooling of drops that strike the wing. 



AERODYNAMIC HEATING 



In viscous fluid flow the losses in mechanical energy 
appear as heat. If all the losses in mechanical energy 
in a fluid element are assumed to appear as heat in tha t 
element, then the total energy of the element will remain 
constant. A surface exposed to "boundary-layer friction, 
"but otherwise thermally isolated, would then "be expected 
to reach a temperature higher than the outside air "by an 
amount determined "by the assumed relation that the excess 
temperature energy of the air stopped at the surf ace equals 
the mechanical energy of the outside air. 

The validity of the preceding 'assumption can be ex- 
amined for cases where the boundary layer is thin. The 
total energy per unit mass of an element of mass of the 

air is — + Jc^T, where u is the velocity in the di- 

rection of flow; J, the mechanical equivalent of heat; 
Cp, the specific heat at constant pressure; and T, the 
absolute temperature. 

It will be shown that the total energy is unchanged 
by changes in pressure and, under certain conditions, by 
viscosity and heat conductivity. If viscosity is neglect- 
ed , 



where p is the pressure and p is- the density of the 
air.. For perfect gases this .expression becomes (if hea 
conductivity is neglected) 



pu du = -dp 




In a thin boundary layer the contribution due to the 
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heat conductivity k is Jk ~ — §• and that due to the 



viscosity |j, is \x -~ ^u |~ where y is the coordi- 

nate perpendicular to the surface. The total energy' con- 
tribution from these sources per unit time is 

dE . d 2 T 3 / du \ , N 

— = Jk + p, — (u ~ ) (1) 

dt dy 2 3y N 3 y / 

When the fluid enters the "boundary layer j these two 
contributions are related if the fluid started from a 

uniform stream because, in this case, — + Jc^T = constant 

2 * 

holds initially, Thus, equation (l) becomes ini t ially 

_ 3 2 / u s \ 9 / 3u\ dE 
_J k ( w ^ f u \ = 

3y 2 \ 2Jc p / dy ^ 3y ; dt 

or . 

\ k /oy 2 V2 / . dt k 

Nov/, if Cp|i,/k, which is called cr or the Prandtl number, 
is unity,' dE/dt is zero and i. u 2 + J Cp- will remain 
constant. 

This analysis is also valid for turbulent boundary 
layers if m- includes the eddy viscosity and k the 
eddy heat conduction. For turbulent boundary layers, 
then, O" will be closer to unity because the eddy vis- 
cosity is equal to the eddy heat conduction divided by 
the specific heat, as in .Reynolds 1 analogy. 



It is readily shown that the distribut ion u* / 2 + 
JcpT = constant is the condition for no heat transfer 
from the surface because in that case 
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A surface exposed to a "boundary layer (in a fluid for 
which cr = l) "but otherwise thermally isolate.d will there- 
fore reach a temperature T 0 + u 0 2 / 2Jc p» ( T o and u o 
are free-stream values,) This result has "been previously 
obtained fo* c ompre s s ible laminar flow along a flat plate 
( reference l) . 

For air, cr is 0.769 (at 0° C) for laminar boundary 

layers and, in this case, — — i s positive; hence dE/dt 

9y 3 2 

is negative and a rise in temperature of the stopped fluid 
somewhat less than the adiabatic rise (u 0 2 /2Jcp) would be 

expected. Fohlhausen (see reference 2, p. 627) calculated 
the temperature rise for the case <T £ 1 and obtained ap~ 

proximately <r 2 (u 0 /2JCp) ; and for (J = 0.769, his re- 
sults give 88 percent of the adiabatic rise. For the tur- 
bulent boundary layer, cr is closer to unity and a some- 
what larger temperature rise should be expected. 

The ITACA 24-inch high-speed wind tunnel (reference 3) 
offers an excellent opportunity to test these theoretical 
conclusions for two reasons: First, high speeds result in 
large, easily measured temperature differences; and second, 
the test-sect ion air has been expanded from atmospheric 
conditions. Thus, the determination of the difference be--, 
tweon surface temperature and atmospheric temperature will 
give a direct measure of the difference between the sur- 
face temperature and .the temperature the surface would 
have had if it had experienced the adiabatic temperature 
rise. This difference being much smaller than the actual 
temperature rise, a considerable gain in accuracy is ef- 
fected by this method. 

Tests were therefore conducted in the NACA 24-inch 
tunnel of an. airfoil prepared as shown in figure 1. Plat- 
inum wires 1 foot long having a diameter of 0.002 inch 
were embedded in balsa inserts in a 5-inch-chord aluminums- 
alloy ITACA 0012T airfoil. The balsa inserts were painted 
and the paint was rubbed to a smooth surface exposing the 
wires. The wires were so placed that one would be in the 
laminar boundary-layer region (at 10 percent of the chord) 
and one in the turbulent layer (at ,70 percent of the 
chord) . 

The surface temperatures were determined by measuring 



NACA Technical Note No, 779 



5 



the resistance of the vires, which w.ere calibrated after 
the wing was tested. .Atmospheric temperatures were 
found "before and after each test polnt : was taken by de- 
termining the wire resistances at a very low tunnel speed 
(30 to 40 mph) . 

In order to check the method of measurement, a wire 
was inserted at the stagnation point of an NACA 0018 air- 
foil and the temperatures were measured over the range of 
speeds. As was anticipated, no significant difference in 
temperature "between the stagnation point and the atmosphere 
was found to exist. 

The results for the NAG A 0012T airfoil at 0° angle 
of attack are plotted in figure 2. The temperature rises 
are nearly equal to the adiabatic, as was expected. The 
fact that the wire in the laminar flow gave a considerably 
lower temperature than the one in the turbulent flow is 
chiefly due to the fact that the pressure at 0.10c was 
lower than at 0.70c. This lower pressure produces a 
greater temperature drop, in the fluid just outside the 
boundary layer and therefore, even if the same fraction 
of this drop were recovered across the boundary layer, a 
lower temperature would be expected, at the lower pressure 
point.. 

In order to obtain a comparison of the test results 
with theory, the temperature rises across the boundary 
layer were computed. The local pressures were obtained 
by calculating the low-speed pressure. distribution by 
Theodorsen's method (reference 4) and correcting for com- 
press ibili ty by multiplying these pressures by 
1/ v'l M 3 ./ Here M , the Mach number, is the ratio of 
the velocity in the test section to the velocity of sound 
in the test section. The local Mach number and the temper- 
ature of the air outside the boundary layer were then cal- 
culated from this pressure distribution for- the laminar 
and the turbulent stations. The- percentage of the local 
temperature drop' outside the bdundary ' lay.er that was re- 
covered in the boundary lay er . i s " shown, in figure 3. It 
should be pointed put that the -points .'pl-otted for Local 
Mach number of the order of 1.0 or '.greater -have little 
meaning because the preceding method of calculating the 
local pressure is invalid in these cases* 

The temperature rises across the laminar boundary 
layer are only slightly smaller than those across the tur- 
bulent layer. The scatter of the points in the turbu- 
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lent layer was also greater than in the laminar layer. 
It was thought that perhaps the boundary 1 layer at the 70- 
percent station might he laminar. part of the time.. This 
possibility was eliminated by noting that a strip -of Ho. 
80 emery placed ' just back of the wire at 0.10c had no 
effect on the temperature rise of the wire at the 0.70c 
s tat ion. 

The temperature rises for high speeds are large 
enough to have an important effect on aircraft icing. In 
particular, airplanes traveling at speeds greater than 
300 miles per hour would probably not experience the most 
severe icing conditions, which occur at atmospheric tem- 
peratures only a few degrees below freezing.. Aerodynamic 
heating is also probably the reason for the. lack of any 
ice formation on the outboard regions of propellers ex- 
cept under conditions of extreme cold and low tip speed. 

If the ice formed and was slung off by centrifugal 
force, the* shear strength of ice indicates that pieces 
less than one-eighth inch thick, which, ought to be light 
enough to stick, should bo. found. 

This phenomenon ought to make it easier to. supply 
enough heat to prevent ice formation on high-speed air- 
planes because a considerable part of the necessary heat 
will be supplied by the boundary-layer friction. 

PATHS' OF WATER DROPS 



The most dangerous ice formations are produced by 
the striking and the freezing of the supercooled water 
drop on the aircraft. Inasmuch as supercooling of these 
drops can remove an appreciable amount of heat, the frac- 
tion of the swept area that is cleared of drops by a wing 
should be of interest in the prevention of ice formation. 
Also, for cases where the drops freeze on impact, as in 
the formation of rime, a relation between the iced area 
and the maximum drop size can be obtained from a knowledge 
of the paths. In order to provide qualitative information 
on these and other points, numerical " calculations of the 
drop paths in air flowing at 200 miles per hour around a 
circular cylinder 1 foot in diameter were made. 

The equation of motion' of a drop is. 



1TACA Technical iTote "!To. 779 



7 



m 



dV 
dt 



(6) 



v/here m is the mass of the drop; V, its velocity; 
D' f the drag of the drop'in an air stream of velocity 
(V - V a ) , v/here V a is the local air velocity and 
(V - V a ) x is a unit vector. The drag, can "be taken from 
data for spheres given "by Zahia (reference 5) for a drop 
of diameter d as 



d = / 28 rn^. 

^ L v 



V a )d 



-0.85 



J— ( 

+ 0.48} i \ p (V - V a ) S U£ 



(7) 



The local air velocity around the forward part of a cyl- 
inder can "be found from the potential function (refer- 
ence 6, p. 30) 



r 2 y 



where XJ is the stream velocity; r, the cylinder radi- 
us; and the coordinates are measured from the center of 
the cylinder. 



The velocities parallel and perpendicular to thxs 
stream direction are 



u = 



chjjr 
9y 



2r 2 y 



2 V 3 



x 8 + y 2 



(x 2 4- y 2 ) 2 



- .M = U 

3x (x2 + ya) 3 



(8) 



The equation of motion of the drop was then numeri- 
cally integrated. A drop of a given size was assumed to 
"be at a given ordinate y 0 and at an abscissa x Q large 

enough that the initial acceleration was small; the drop 
was also assumed to have the free-stream velocity. The 
local air velocity V a was then obtained from equation 
(8) and the air force and the acceleration from equations 
(7) and (6). The drop was assumed to move with this ve- 
locity and acceleration for a short time interval At 
and new coordinates and a new velocity v/erc computed after 
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this interval. The process v/as repeated until the drop 
either hit or clearly missed the cylinder. In this way, 
the maximum initial ordinate that a drop could have and 
still strike the cylinder v/as found. Thus, the area 
swept clear of drops v/as. found and is plotted in figure 4» 
From the same calculations, the region of the cylinder 
struck "by the drops v/as obtained and is plotted in figure 
5. This, region would "be the area covered with ice if the 
drops were of uniform size and froze on contact. 



Langley Memorial Aeronautical Laboratory-, 

National Advisory Committee for Aeronautics., 
Langley Field, Va. , July 24, 1940. 
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Figure 1.- Cross section of aluminum- alloy NACA 0012T airfoil 
used in tests,. 
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Figs. 3, 3 
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Figure 2.- Variation of observed temperature rise with air speed. 
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Figure 3.- Percentage of adiabatlc rise recovered in the boundry layer. 
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Figure 4.- Variation of swept area cleared of drops with drop diameter for a right circular 
cylinder. 
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Figure 5.- Variation of area of right circular cylinder struck by drops with drop diameter. 



